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We theoretically analyzed inelastic effects in the electron transport through molecular junctions 
originating from electron-vibron interactions. The molecular bridge was simulated by a periodic 
chain of identical interacting hydrogen-like atoms providing a set of energy states for the electron 
tunneling. To avoid difficulties inevitably arising when advanced computational techniques are 
employed to study inelastic electron transport through multilevel bridges, we propose and develop a 
semiphenomenological approach. The latter is based on the Buttiker's dephasing model within the 
scattering matrix formalism. The advantage of the proposed approach is that it allows to analytically 
study various inelastic effects on the conduction through molecular junctions including multilevel 
bridges. Here, we apply this approach to describe features associated with electron energy transfer 
to vibrational phonons which appear in the inelastic tunneling spectra of electrons. In the particular 
case of a single level bridge our results agree with those obtained by self-consitent calculations carried 
out within the nonequilibrium Green's functions method validating the usefulness of the suggested 
approach. 

PACS numbers: 31.15.at, 72.10.Di 
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I. INTRODUCTION 

Molecular electronics is currently recognized as one of 
the most promising developments in nanoelectronics, and 
a significant progress had been seen in this field dur- 
ing the past decade [J, The key element and basic 
building block of molecular electronic devices is a molec- 
ular junction consisting of two conducting leads linked 
by a molecular bridge. Conduction through molecular 
junctions is being intensively studied both theoretically 
and experimentally, and main transport mechanisms are 
presently understood 0, However, the progress of 
experimental capabilities in the field of molecular elec- 
tronics unceasingly brings new challenges urging further 
development of the theory. 

In practical molecular junctions electron transport is 
always accompanied by nuclear motions in the surround- 
ings, so the conduction is affected by the coupling be- 
tween electronic and vibrational degrees of freedom. To 
analyze vibration-induced effects in the electron trans- 
port through molecular bridges one assumes that molec- 
ular orbitals are coupled to phonons describing the vi- 
brations, and this brings an inelastic contribution to the 
electron current. Over the past decade theoretical stud- 
ies of inelastic vibration-induced effects in the electron 
transport through molecules and other similar nanosys- 
tems (e.g. carbon nanotubes) were carried out by a large 
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number of authors (see Refs. 0, H, H, 0, i, i, [13, 
12,[l3|)- Also, manifestations of the electron-vibron in- 
teractions were observed in the experiments using vari- 
ous molecules, carbon nanotubes and fullerenes as linkers 
in the junctions [H [H [H, [13, [11, [3 S HH, [12, El] . 



Vibration-induced features were viewed in the differential 
molecular conductance dl/dV (J, V respectively being 
the electron current through the junction, and the bias 
voltage applied across the latter) and in the inelastic tun- 
neling spectra of electrons (Pl/dV^. Sometimes even the 
I — V curves themselves exhibit an extra step originating 
from the electron-vibron interactions, as was reported in 
the Ref. 17[ where the electron transport through a sin- 



gle hydrogen molecule placed in between platinum leads 
was explored. 

Here, we focus on theoretical studies of the inelastic 
electron tunneling spectra (lETS) which is the rich and 
diverse problem attracting an unceasing interest of the 
research community for a long while 0, [o, S 13 23, 2^ ■ 
A systematic framework for such studies is based on 
the nonequilibrium Green's functions formalism (NEGF). 
Within the NEGF the electron transport through a two 
terminal junction including a linking molecule (presented 
as a set of energy levels Ei ) coupled to a set of phonons 
with the frequencies f2i could be properly analyzed. 
However, the extreme complexity of the general equa- 
tions for the relevant interdependent electron and phonon 
Green's functions limits the practical applicability of 
NEGF to simple models. Currently, specific NEGF based 
results for lETS are mostly obtained for a model where 
the molecule is simulated by a single orbital (a single- 
state bridge) which is coupled to a single vibrational 
mode with the frequency Adopting this simple model, 
it was shown that a vibration-induced signal appears in 
the lETS near the phonon excitation threshold eV = hfl. 
The signal looks as a peak, dip or derivative-like feature 
in the voltage dependence of d^I/dV^. Its specific shape 
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is sensitive to the junction characteristics such as the 
electronic state energy, electron-phonon and molecule-to 
leads coupling strengths and the vibronic frequency (see 
the recent review and references therein [7|). 

In the present work we analyze the lETS using more 
realistic model to mimic a molecule linking the terminals 
in the junction. We simulate the molecule by a peri- 
odical chain of identical hydrogen-like atoms with the 
nearest neighbors interaction. The first and the last sites 
in chain are coupled to the leads. Such a model was 
fist proposed by D'Amato and Pastawski [2^, and af- 
terwards it was repeatedly used to study electron con- 
duction through molecules. The chosen model gives us 
to describe a multisite bridge for the electron transport. 
The number of electronic energy states on the bridge 
equals to the number of the bridge sites [l^], and the 
states are situated within the energy range whose width 
is determined by the coupling between the sites. To 
avoid computational difficulties arising when the con- 
sistent NEGF is employed, we chose a semiphenomeno- 
logical approach based on the combination of the latter 
with the Buttiker dephasing model [2^. This enables 
us to significantly reduce computational effort. Earlier, 
the similar approach was successfully applied to analyti- 
cally analyze dissipative electron conduction through sin- 
gle site molecular bridges [1^. Within the chosen ap- 
proach we analyze the vibration-induced features in the 
(Pl/dV^ versus V spectrum, and we show that the spec- 
trum may significantly differ from the spectra typical for 
single site molecular bridges. Besides the peak/dip or 
derivative-like feature at the phonon excitation thresh- 
old, extra lETS features may appear, whose location, 
size and shape strongly depend on the number of sites 
and site - to site couplings in the bridge. 



II. MODEL AND MAIN EQUATIONS 

An important advantage of the phcnomcnological 
model for the incoherent/inelastic quantum transport 
proposed by Buttiker ^28^ is that this model could be 
easily adapted to analyze various inelastic effects in the 
electron transport through molecular junctions. Within 
this model one may analytically study such effects avoid- 
ing inevitable technical difficulties which arise when the 
advanced methods are applied to compute characteristics 
of the electron transport through multisite (and multi- 
level) molecular bridges. Following Buttiker's approach, 
inelastic effects are accounted for by means of electron 
reservoirs attached to the bridge sites. Electrons could be 
scattered into these reservoirs where they participate in 
dephasing/inelastic processes before their returning back 
to the bridge. 

In further calculations, we mimic a molecular bridge 
linking two electrodes as a periodic chain of identical sites 
including 2m— 1 elements. We set the site ionization en- 
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FIG. 1: (Color online) Schematic drawing illustrating electron 
transport through the junction including a chain of 2m — 1 
identical sites. The dissipative electron reservoir is attached 
to the middle site. The semicircles represent left (L) and 
right (R) electrodes. 



ergies Ei = Eq and nearest neighbors coupling strengths 
= — /3- The schematic of the adopted model 

is shown in Fig. 1. The first and the last sites in the chain 
are attached to the leads and remain fixed. We suggest 
that the chain vibrates at a fundamental frequency fl, 
and we study the electrons interaction with the corre- 
sponding transverse mode. In this case, the maximum 
displacement of the chain sites occurs in the middle of 
the latter. Therefore we assume, as a first approxima- 
tion, that only the middle site in the chain is coupled to 
the vibrational mode. To describe this coupling within 
the Buttiker approach, we put the middle site in con- 
tact with an electron reservoir. The electrons scattered 
in there could absorb/emit phonons and then reappear 
on the bridge. 

In the figure 1, the squares represent barriers sep- 
arating sites included into the chain and the triangle 
imitates the scatterer coupling the middle site to the 
reservoir. The scatterer is described by the scatter- 
ing matrix s which relates outgoing from the scatterer 
wave amplitudes a'l , 02 , 03 , 04 to the incorning ampli- 
tudes ai, a2, as, a4. The latter has the form [28[: 
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(1) 



Here, the parameter e takes on values within the range 
[0, 1]. This parameter determines the probability for an 
electron to be scattered, and it was introduced in Ref. 
[2^ as a phenomenological characteristic. In further 
analysis, we suggest to express e in terms of relevant 
energies characterizing electron transport through the 
molecular junction (see the next Section). This provides 
means to better describe the transport processes within 
the chosen scattering matrix formalism. 
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To simplify the following calculations of the electron 
transmission, we assume that an electron could be in- 
jected at the system shown in Fig. 1 and/or leave from 
there solely via four channels indicated in the figure. 
Then incoming particle fluxes Jk are related to those 
outgoing from the junction (J 'A by means of the 4x4 
transmission matrix T 12 8l. 1291 : 



fc 



TjkJk'i 



J<k<4, 



(2) 



where Ji, J[ = Jl, J'l correspond to the left side of the 
junction and the fluxes ^2,^2 = JRtJ'r ^^'^ associated 
with its right side. To provide charge conservation in the 
system, the net particle flux in the channels 3,4 must 
equal zero. 

Assuming the electron transport from the left to the 
right, the electron transmission is the ratio of particle flux 
outgoing from the right end of the junction 2 and that 
incoming to its left end 1. To arrive at the expression 
for T{E), one must solve the system of linear equations 
(121). This gives: 



T{E) = T21 



Ki{E)K2{E) 
2 - R{E) 



(3) 



Here, 



(E) - Tai + T41 ; K2 (E) = T23 + T24. (4) 



R{E) = T33 + + T43 + T44, 



(5) 



The matrix elements Tjk introduced in the Eq. Q are 
closely related to the elements of the scattering matrix 
S describing the whole system {Tjk = \Sjk\'^). The lat- 
ter expresses outgoing from the system wave amplitudes 
6^, 6^, 03, 04 in terms of the incident ones 5^^, 5^^, 03, a4. 



The matrix elements Sjk are computed below. To reach 
the middle site in the molecular chain, an injected elec- 
tron must flrst tunnel from the left electrode to the chain, 
and then it tunnels through m — 1 barriers separating 
the sites. The tunneling through this set of barriers may 
be described by the 2x2 matrix Wl whose elements 
are expressed in terms of transmission and reflection am- 
plitudes for the barriers: 



Wr = 



''M 



Knr + Kn) =1- (6) 



Assuming all barriers separating the sites to be identi- 
cal, we may write the following expression for this matrix: 



Wr = 



1 










1 








(7) 



where ,t and r are the transmission and reflection 

coefficients for the barrier between the left electrode and 
the chain, and for a single barrier separating the sites 
on the left half of the chain, respectively. The electron 
tunneling through the right half of the chain could be 
described in a similar way by using the matrix Wr : 
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and reflection {r^ 



\ 
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^ / 

(8) 

) coeffi- 



where the transmission {t 
cients characterize the barrier between the chain and the 
right electrode. 

Combining Eqs. H]), © and ^ we obtain the follow- 
ing expression for the scattering matrix S : 



L 

rn 



a 



,,L ^.R 



(9) 



m m I 



Here, Z — u^u^ — a^v^v^] a = \/T~~e; (3 ~ -^e. For a 
single site bridge symmetrically coupled to the electrodes 
(to = 1, =t^ = t'] =r^ = r') Um = Vm = 

—r'/t\ and the Eq. ^ agrees with the corresponding 
result of Ref. Inserting the expressions for matrix 

elements Tjk determined by Eq. Q into the Eq. ([2]), 
we may present the electron transmission in the form: 



T(E) = ^ ^ ^ '^ (10) 

2K+a2(^2^_l)]2 

where Wm could be expressed in terms of matrix ele- 
ments of the matrices Wl,r- 

Our next step is to find the relation of the quantities 
Wm to the relevant Green's functions. Disregarding for 
a while the effects of nuclear motions, we may write the 
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following effective Hamiltonian for the molecular bridge: 



(11) 



Here, the first two terms describe the bridge itself. 
Within the adopted approximation their nonzero matrix 
elements between states and \l) corresponding to 
the kth and Ith sites on the molecular chain have the 
form (ffo)fci = EqSm and {Hi)ki = + 
The remaining terms represent self-energy corrections 
arising due to the coupling of the leads to the chain. 
Supposing that only the first and the last sites of the 
molecular bridge are directly coupled to the correspond- 
ing leads, we can write 

{HL)kl ~ SklSkl^L', {HR)kl = SklSk,2m-l^B.- (12) 

The self-energy terms could be written in the form: 



LM\2 



kr 



E — er + i(T 



(13) 



In the expression (jl3p summation is carried out over the 
states on the left/right lead, is the single electron 
energy for electron state r in the lead, r characterizes 
the coupling strength between the lead and the molecular 
chain, and tr is a positive infinitesimal parameter. 

The electron transmission for the elastic coherent 
transport through the considered junction is presented 
as: 



T = rLTR\Gi, 



(14) 



where Tlh = — 2Im(I]L^/j) , and the relevant Green's 
function matrix element is given by: 



Gi,2m-i = {1\E-Heff\2m-1). 



(15) 



Further, we proceed within the wide band approxima- 
tion for the leads. Within this approximation self-energy 
terms (|13|) do not depend on the electron tunnel energy 
E. Also, we disregard the real parts of Sl,^. Then the 
expression for the Green's function may be obtained in 
the form described in the Ref. 27|, namely: 



Gl.,2m-l{E) 



(A + c)2('"-i) (A + c + «r)2 - (A - c)2("-i) (A - c + «r)2 ' 

(16) 

Here, we introduced denotations: 
A = Eo-E; C = VA2-4/32; r = Tl+Tr. (17) 



In the absence of phonons {a = 1), the expression ([10] 
for the electron transmission takes on the form: 



T{E) = 



1 



(18) 



where 9f-^ — l/wf^^. Comparing this result with Eq. ([M]) 
we obtain: 



1 



9n 



(19) 



where the function gm{E) is closely related to the 
Green's function of the chain, namely: 



9r 



,{E) = ^TLTR\Gl,2m~l{E)\ 



(20) 



Using the result of ^1% we may rewrite the general ex- 
pression for the electron transmission (fTO|l as follows: 



T{E) = 



g„,(^)(l + a") [1 -f g,r,{E) -a\\~ g„,(^))] 
[l + g^{E)+a^{l-gm{E))\'' 



(21) 

This expression is used as a starting point in the further 
analysis. 



III. RESULTS AND DISCUSSION 

In Fig. 2, we show the energy dependencies of T{E) 
calculated assuming the parameter e to be independent 
of the energy E. As expected, instead of a single maxi- 
mum at E — Eq we get a set of peaks located within the 
interval Eq — 2/3 < E < Eq + 2/3. The number of peaks 
equals to the number of sites in the chain. These peaks 
correspond to the energy levels of the molecular bridge 
within the adopted model. The coupling of the chain 
ends to the leads affects the minimum values of transmis- 
sion, especially near E = Eq. As the coupling strength- 
ens, the transmission minimum values increase. When 
the reservoir is attached to the chain (e > 0), it causes 
a decrease in the range of the transmission variations. 
However even within the strong dephasing/inelastic limit 
(e = 1) the transmission peaks remain well distinguish- 
able for multisite molecular bridges. On the contrary, the 
peak in the transmission through the single site bridge 
a^ears completely washed out under strong dephasing 

In the considered case of the vibration-induced inelas- 
tic contribution to the molecular conductance, the pa- 
rameter e must depend on the quantity 'Sph{E), which 
is the self-energy correction arising due to the electron- 
phonon interaction. Comparing the results for the elec- 
tron transmission obtained basing on the phenomenolog- 
ical Buttiker's model with those derived within NEGF 
formalism, it was shown that for a single state bridge e 
is rather simply related to Sp/j, namely (29j : 



ph 



r + r 



(22) 



where Tph — — 2Im(Ep^), and F describes the coupling 
of the bridge to the leads. Presently, we are considering a 
multisite bridge but we assume that only one site in the 
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FIG. 2: (Color online) Energy dependencies of the electron 
transmission through the junction (top and middle panel) and 
the parameter e (bottom panel). The curves are plotted 
assuming e = (dash-dotted lines), e = 0.5 (dashed lines), 
and e = 1 (solid lines) for N = 5 (top panel) and N = 11 
(middle panel). The curves in the bottom panel are plotted 
at iV = 11 (dash-dotted line), = 9 (dashed line), N^7 
(solid line) for M = O.SeV, HQ = 0.22eV. In plotting all 
curves included in the figure it is assumed that r_t ~ Tr = 
0.2eV, P = 0.6el/, Eo = 0. 

middle of the latter is coupled to the vibrational mode. 
Therefore we approximate e by the expression similar to 
Eq. (|22p where the term T in the denominator is replaced 
by the energy A , which depends on both the bridge ends 
coupling to the leads (F) and site to site couplings on the 
bridge itself (/3) . For sufficiently long molecular bridges 
[N > 10) one may expect the value of A to be mostly 
determined by the coupling between the sites. 



The form of Tph varies depending on the characteris- 
tics of the electron-phonon interaction in the system, and 
on the kind of phonons contributing to these interactions. 
For instance, the electron coupling to the thermal phonon 
bath results in the expression for the Tph , which strongly 
differs from that obtained while considering the electron 
coupling to a few vibrational phonon modes. Here, we 
assume that a single vibrational mode with the frequency 
n exists in the system and that the temperature is low 
{kT ^ hil.,T,(3). Then Tph could be presented in the 
form Q: 

Tph{E) ^ttM^ I J (hjpph{oj)pei{E + hui) 

+ J dujpph{uj)peliE - lThj)> (23) 

where M is the electron-phonon coupling strength, pl.r 
are the chemical potentials for the leads, and pR < E < 
PL. The phonon density of states pph{ijj) included in 
this expression reveals a sharp maximum at lo = Q, : 

Pph {^)^\ 7 1^^ 2 ■ (24) 

Here, the linewidth is of the order of the thermal 
energy kT. The electron density of states (DOS) of the 
middle site of the chain is described by the expression 

Pei{E) = ln\[Gmm{E)]. Within the adopted model 

TT 

we obtain: 

Gmm{E) — — 

(A + O^'HA + c + ^r) - (A - c)"'~^(A - c + »r) 
(A + O^'KA + c + *r) + (A - C)™-HA - C + «r) ■ 

(25) 

Now, we may apply the obtained approximation for the 
parameter e to compute the vibration-induced electron 
transmission, and then we employ the latter to calculate 
the electron current through the junction. We accept the 
usual Landauer expression for the current: 

^ = ^ / nmUE) - fR{E)]dE (26) 

where fL,R{E) are Fermi distribution functions for the 
leads with chemical potentials pl,b.- For a symmetri- 
cally coupled junction pl,b. = Ep ± ■^V, where Ep is 
an equilibrium Fermi energy of the junction. The cur- 
rent given by Eq. (|26p includes the vibration-induced 
inelastic contribution. The latter is brought in there by 
means of the expression ()21|) for the electron transmis- 
sion T{E). As shown above, the parameter = 1 — e in 
the Eq. (PT|) depends on the electron-vibron interaction, 
and this leads to the occurrence of the inelastic term in 
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the current along with the elastic one. The proposed for- 
malism enables us to qualitatively analyze the effects of 
the electron- vibron interactions in the lETS for multisite 
molecular bridges. 

The electron- vibron interaction caused self-energy Tph 
strongly depends on the tunnel energy E, and this re- 
sults in a well pronounced energy dependence of the pa- 
rameter e shown in the bottom panel i Fig. 2. We re- 
mark that e remains zero for energies within the range 
Eq — MI < E < Eq+HH. This fits into the accepted model 
where electron- vibron interactions with a single phonon 
mode are assumed to mainly cause inelastic effects in the 
electron transport. An electron on the bridge may vir- 
tually absorb and emit phonons of the energy hfl thus 
creating metastable states around the energy states of the 
molecular chain. In general, the bridge energy states are 
shifted from their original positions due to the electron- 
phonon interactions but in the present consideration we 
disregard these shifts for simplicity. Within our model 
for the bridge one of the bridge states has the energy 
Eq, the same as the ionization energy of noninteract- 
ing sites. Being initially at this energy level, an electron 
could make transitions to the states E — EqzL nhil as a 
result of its interaction with the vibrational mode. So, 
change in the energy of such electron originating from the 
interactions with vibrational phonons cannot be less than 
hil. Otherwise, the electron transport must be coherent 
and elastic (e = 0). At E ~ Eq ± hft, the parameter 
e sharply increases revealing peaks, which correspond 
to the appearance of the metastable vibration-induced 
states around E — Eq. The same consideration could 
be applied to another electronic states on the bridge. As 
expected, e versus E plot displays a sequence of peaks 
whose number is related to the number of elements in- 
cluded in the chain. All these features are located in the 
energy range Eq — 2j3 < E < Eq + 2/3, and e again 
becomes zero outside of this range, as shown in the Fig. 
2. 

Due to the interaction with the phonon mode, the elec- 
tron transmission gets extra features which are revealed 
in the observable characteristics of the electron current 
through molecular junctions. The specifics of these man- 
ifestations is determined by the relative values of the rele- 
vant energies. At low temperatures {kT <C ^l,r, (3, M) , 
the relation of the bridge-to leads coupling strengths 
^L.R and electron- vibron coupling constant M is espe- 
cially important. At weak coupling of the bridge to the 
leads, the electron-vibron interaction splits every single 
peak in the coherent electron transmission (located, for 
instance, at E = E) into a set of smaller peaks asso- 
ciated with vibrational levels at E = E + nhil which 
could be resolved when hfl < F^ j^,/?. Phonon-induced 
peaks in the transmission could give rise to the steps in 
the I — V curves and rather sharp features in the dif- 
ferential conductance and lETS plotted against the bias 
voltage applied to the junction. This was first theoreti- 
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FIG. 3: (Color online) Inelastic electron spectrum d^I/dV 
against V for a molecular junction with a single site (top) and 
multisite (bottom) bridge. The curves are plotted assuming 
Fl = Ffl = 0.2eV, P = QAeV, Eo = 0, hQ. = 0.22eV, M = 
0.3eV, N = 7 (bottom line). 



cally shown by Wingreen et al '3u| and then confirmed 
in some other works (see Refs. ^7,, ^]) for the case of a 
single state bridge. Here, we concentrate on another sit- 
uation, assuming that the bridge coupling to the leads is 
rather strong: F^^ij ~ M. In this case inelastic effects do 
not cause significant changes in the I—V characteristics. 
However they cause the appearance of a distinct signal in 
the inelastic tunneling spectrum at eV ~ HQq similar to 
that occurring on a single-site bridge interacting with the 
vibrational mode. The shape and magnitude of the sig- 
nal are sensitive to the characteristics of the considered 
system. 

We remark that our semiphenological approach brings 
the result which agrees with those earlier obtained in sev- 
eral theoretical works basing on the consistent NEGF for 
single state molecular bridges. The size of the deriva- 
tive -like lETS signal at eV « H^l depends on the num- 
ber of sites in the bridge chain as shown in the figures 
3 and 4 (top left panel). For iV = 1 (a single-site 
bridge), the signal size takes on value of the same order 
as that reported by Galperin et al 0] for a symmetrically 
coupled junction with relatively close values of the rele- 
vant parameters. For multisite bridges the lETS reveals 
more signals which indicate inelastic electron transport 
through the channels associated with the bridge energy 
levels situated apart from E = Eq. An example of the 
vibron-induced lETS for a multisite bridge {N — 7) is 
presented in the right panel of Fig. 3. Some lETS signals 
are separately displayed in Fig. 4 to show their shape and 
size for different values of N. 

Our analysis shows that the signal indicating the 
phonon excitation threshold is accompanied by another 
feature appearing at eV ~ 2MI, and this happens for 
multisite and single-site bridges in the same way. The 
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FIG. 4: (Color online) The lETS features in (fl/dV'^ for 
the several resonant levels model (multisite bridge) plotted 
assuming Tl=Tr = 0.2eV, 13 = 0.4eV", hQ. = Q.22eV, M = 
0.3eV, £o = 0. Top panels: iV = 11 (dash-dotted lines), 
A'^ = 3 (dashed lines), A'^ = 1 (solid lines). Bottom panels: 
iV = 11 (dash-dotted lines), N = 7 (dashed lines), = 3 
(solid lines). 



multiplication of lETS signals is known for molecular 
junctions with weak coupling of the bridge to the leads, 
and it shows the contribution of higher phonon harmonics 
into the electron transport. However, this explanation is 
hardly justified for the presently considered junction with 
rather strong coupling of the bridge to the leads. To clar- 
ify the nature of the obtained lETS signal duplication, we 
computed (Pl/dV^ for a junction with a slightly asym- 
metric coupling. It appears that the extra signal position 
varies depending on the ratio Fi/Ffl, (see Fig. 5), and 
the higher is the asymmetry, the closer it moves to the 
phonon excitation threshold. In general, the extra sig- 
nal appears at (1 -I- TL/TR)eV ~ Ml. One may expect 
that within the scanning tunneling microscopy (STM) 
junction configuration (Fl <C T b) this signal would be 
superimposed on that indicating the excitation thresh- 
old for the vibrational phonons. This gives grounds to 
conjecture that the effect of the electron- vibron interac- 
tions on the lETS is twofold. The very opening of the 
channel for inelastic transport gives rise to the signal in 
the inelastic tunneling spectrum of electrons. This sig- 
nal always occurs at the phonon excitation threshold re- 
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FIG. 5: (Color online) The effect of asymmetric coupling of 
the bridge to the leads on the lETS feature in (fl/dV'^. The 
curves are plotted for N = 7, Tl = Tr = 0.2eV (solid 
line) and Fl = O.lSeV, Vr = 0.22eV (dashed line). The 
remaining parameters take on the same values as in the Fig. 
4. 



gardless of the characteristics of the bridge coupling to 
the leads. Another signal indicates that the correspond- 
ing metastable electronic state appears in the conduc- 
tion window between /i^ and /i^j , and its position is 
determined by the bias voltage distribution in the junc- 
tion. The two coincide within the STM configuration, 
and their maximum separation happens in symmetrically 
coupled junctions and equals Kl. Similar signal duplica- 
tion arises due to the opening of inelastic transport chan- 
nels associated with every electronic state of the bridge 
chain. 

Finally, we remark that the adopted computational 
approach proved itself workable. It allows us to repro- 
duce the principal vibration-induced inelastic effects in 
the electron transport through molecular junctions and 
to save computational efforts and time. In the partic- 
ular case of a single-site bridge linking the leads in the 
junction, the proposed approach gives results which agree 
with those obtained using advanced computational for- 
malisms. At the same time, the adopted model enables 
us to analyze inelastic electron transport through mul- 
tisite bridges simulating practical molecules. One may 
generalize the proposed computational method to include 
into consideration a set of vibrons with different energies. 
For this purpose one must include into the model extra 
dephasing reservoirs attached to the bridge sites. Also, 
electron structure of a practical molecule may be used 
in calculations of the bridge Green's functions, and this 
could help to further improve the results, and make them 
more useful and convenient for comparison with experi- 
ments. 
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